Bean-Livingston barrier enhancement on nodal surface of the d-wave superconductor 
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Vortex entry into (110) oriented YBa2Cu307_a; films has been studied by tunneling into Andreev 
- Saint- James bound states, whose energy is shifted by surface currents. At low temperatures, the 
characteristic field for vortex entry has been found to increase up to values several times higher than 
that of the Bean-Livingston entry field for conventional superconductors, in agreement with recent 
theoretical predictions. 

PACS numbers: 74.25.Ha, 74.50.+r, 74.60.Ge, 74.72.Bk, 74.76.Bi 
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As shown by C. Bean and J. Livingston [l[ vortex en- 
try in a Type II superconductor submitted to a magnetic 
field parallel to its surface can be delayed beyond the field 
H c i where the mixed state becomes thermodynamically 
stable in the bulk. Vortex entry may not occur up to a 
field of the order of the thermodynamical critical field H c 
where the free energies of the normal and full Meissner 
states become equal. The delayed vortex entry comes 
from the attraction between a vortex and its anti-vortex 
image. The vortex is however also submitted to a repul- 
sive force from Meissner surface screening currents (Fig. 
[TJ At low applied fields these Meissner currents are weak, 
the vortex - anti- vortex attraction wins and a surface en- 
ergy barrier prevents vortex penetration in the bulk. At 
high fields, Meissner currents are strong, the repulsive 
interaction wins, the surface barrier disappears and vor- 
tices can penetrate freely. The attractive and repulsive 
interactions balance each other at a field of the order of 
H c . In practice, vortex penetration occurs in general at 
a field lower than H c , as the sample surface roughness 
can locally enhance the surface field value above that of 
the applied field because of demagnetization effects. 

Recently, Iniotakis et.al. [2j have predicted than in 
a d x 2_ y 2 wave superconductor having its surface normal 
parallel to a nodal direction, vortex penetration for fields 
applied along the z direction (the c-axis perpendicular 
to the CuO planes in the cuprates, see Fig. [2]) can be 
delayed well beyond the field of first vortex entry H s . 
The increase of the barrier, being particularly significant 
at low temperatures, is due to is due to strong surface 
currents carried by Andreev - Saint- James bound states 
3- This current is paramagnetic and flows in the op- 
posite direction with regards to that of the diamagnetic 
screening Meissner current. A schematic illustration of 
the currents is presented in Fig. Q] One must notice that 
the Fogelstrom paramagnetic current flows at a distance 
of the order of the coherence length £ from the surface 
while the Meissner screening current flows at the scale 
of the penetration length A from the surface. The effect 
of the paramagnetic current on vortex entry is to reduce 
the repulsive force driving the vortex into the bulk and 




FIG. 1: Schematic representation of a vortex and its anti- 
vortex image (dotted) together with the sum of Meissner dia- 
magnetic screening current and the Fogelstrom surface para- 
magnetic current. The vortex anti- vortex interaction is at- 
tractive, the effect of the diamagnetic Meissner current is re- 
pulsive and that of the Fogelstrom paramagnetic current is 
attractive. Dashed and full arrows indicate the directions of 
the paramagnetic and diamagnetic currents. At low tempera- 
tures, the sum of both currents becomes paramagnetic on the 
scale of the coherence length £. 



away from the surface. 

In this work, we present a study of the temperature- 
magnetic field dependence of the Bean - Livingston sur- 
face barrier for vortex entry on the nodal surface of 
slightly overdoped YBa 2 Cu 3 7 _ :l; (YBCO) thin films. 
Our results basically confirm the theoretical predictions. 

Thin films of YBa 2 Cu307_ :! ; (YBCO) were grown by 
DC off-axis sputtering on substrates of SrTi03 with (110) 
orientation. A buffer layer of PrBa2Cu307_ :r , was de- 
posited using RF off-axis sputtering on top of the sub- 
strate [J] in order to reduce the (103) orientated grains. 
The YBCO growth time was 3 hours corresponding to 
a thickness of 2400A. Normal metal - Insulator - Super- 
conductor (NIS) junctions were prepared by placingthin 
Indium rods on top of the surface of the thin filmja, [|| . 
Oxidation of the Indium to InO at the contact area gives 
the Insulating layer of the NIS junction. Figure [2] shows 
the NIS junction orientations with regards to the rf-wave 
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FIG. 2: Schematic illustration of the experimental setup. The 
CUO2 planes with the rf-wave order-parameter are scathed. 
The thin film surface is [110] oriented. The magnetic field 
is applied parallel to both the surface and the c-axis. The 
round rod represents the planar tunnel junction (Indium in 
grey, InO in black). 



order parameter. The samples discussed in this work 
were all slightly overdoped with a T c downset of 89K. 

Andreev - Saint-James zero energy surface bound 
states (ASJ bound states) are known to develop at sur- 
faces perpendicular to a node direction, such as (110) 
surfaces in YBa2Cu307_a; @, @|- In a tunneling experi- 
ment, these states show up as a Zero Bias Conductance 
Peak (ZBCP), which splits due to a Doppler effect under 
superfluid currents flowing parallel to the surface 0] as 
occurs in the geometry show in Fig. [3] We have used 
this split 25 to investigate the field and temperature de- 
pendence of the barrier against vortex entry. 

In films thicker than the London penetration length 
as is the case for the samples investigated here, the field 
evolution is known to be hysteretic, being larger in in- 
creasing than in decreasing fields [!, 0, 0|. The split 
in increasing fields has been interpreted as due to the 
Doppler effect induced by Meissner currents 0, [ll[ while 
that in decreasing fields has been shown to be a field 
rather than a current effect [l(|. A third case is that 
of a field cooled condition, in which there should be no 
Bean-Livingston effect. A set of three such measurement 
is shown in Fig. [3] In order to cancel out the field ef- 
fect, we have characterized currents resulting from the 
Bean-Livingston barrier by the difference between values 
of S measured in increasing field after cooling the sam- 
ple in zero field and field cooled conditions, (S up — Sf.c.)- 
This quantity has been measured as a function of field 
and annealing temperature T a as described in the next 
paragraph. 

In the absence of any Bean-Livingston barrier effect, 
(Sup — Sf.c.) should reach a maximum value at H = 
H Cl . In the presence of a barrier against vortex entry, 
(S U p — Sf.c.) should reach its maximum value at the field 
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FIG. 3: dl/dV at 2.5 Tesla for increasing, decreasing and field 
cooled conditions. Defining 28 to be the distance between the 
two peaks, we note 5 is maximum for increasing field, min- 
imum for decreasing field and has an intermediate value for 
the field cooled measurement. Inset shows the measurement's 
full scale (±20mV) for the decreasing field condition. The 
YBa2Cu307_a; coherence peaks are located at ~ l6mV. 



Hg where the barrier disappears (following the notations 
by C. Bean and J. Livingston [l[). Tunneling measure- 
ment at relative high temperatures have thermal smear- 
ing that makes it impossible to follow the ZBCP split- 
ting. In order to follow the temperature dependance of 
the barrier by tunneling measurements and to eliminate 
the effect of thermal smearing on the tunneling charac- 
teristic, we have pursue the following procedure in deter- 
mining the temperature dependence H S (T). After cool- 
ing in zero field down to 1.4K, the field is raised up to a 
value H*. The temperature is then raised up to a value 
T a (hence annealing temperature) for a time sufficient to 
reach equilibrium (2 minutes), then cooled down again 
to 1.4K for a measurement of S up (H* ,T a ). In order to 
return to the virgin conditions (zero field cooled from 
above T c ), the temperature is then raised up to 105K, 
the field is lowered to zero, and a new measurement is 
carried out at a different value of T a . Once several values 
of T a have been explored, a new cycle is started with a 
different value of H*. In order to check if the performed 
field-thermal cycle damaged the sample or the junction, 
the tunneling characteristic was measured every time a 
virgin state was reached and also at the increasing field 
state at the previous H* value and compared to the tun- 
neling characteristics obtained at the beginning of the 
measurements. No change in the tunneling characteris- 
tic occurred during the measurements. 

Fig. |4] summarizes our results for (S up (T a ) — Sf.c.) a s 
a function of H * at different annealing temperatures. At 
T a = 20K and higher, the curves present a broad max- 
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FIG. 4: The hysteresis in the ZBCP splitting value {5 U p{T a ) - 
S F .c.) for T a = 1AK,20K,40K and 50K. At 1.4K the hys- 
teresis increases until the splitting smears into the background 
at fields higher then 3 Tesla. When annealing to higher tem- 
peratures, a broad maximum for the hysteresis appears. For 
T a = 50-ftT, the broad maximum is already achieved at field of 
about 1 Tesla. Solid arrows are guide to the eye. 



imum at a field Hm that goes down as the annealing 
temperature goes up, as expected. At the lowest temper- 
atures, no clear maximum can be seen up to the highest 
field where the value of 6 up (T a ) can be identified. 

Ideally, the field dependence of (5 up (T a )~ Sf.c.) should 
show a sharper maximum than that which we observe 
since a reduction in surface currents and hence of S up is 
expected once vortices start to penetrate. Yet, we ten- 
tatively identify the field at maximum with Hs- What 
is clear from the data is that surface currents keep in- 
creasing up to a field that is much higher than Hc 1 (of 
the order of 1 • 10~ 2 Tesla) at all temperatures, and than 
He (of the order of 1 Tesla [2, [12] ) at low tempera- 
tures, this basically confirms the prediction of Iniotakis 
et.al. In order to check whether thermal fluctua- 

tions are the reason for the continuous change in 5 up , a 
set of measurements was done with much longer time of 
temperature annealing - up to 6 hours. No difference be- 
tween the short and long annealing time was found and 
hence - the thermal fluctuation is ruled out from being a 
reason for this behavior. 

Another way to represent the results is shown Fig. [5l 
where the dependence of of S up (Ta) — 5 up (lAK) is dis- 
played for a fixed field H* as a function of ^f. S(H*) 
is reduced progressively from its value in increasing field 
down to its field cooled value as a function of the an- 
nealing temperature ([5]). This reduction can be fitted to 
an exponential decay, giving a temperature scale T S (H*) 
that decreases as H* is increased. Instead of this progres- 
sive reduction, one would have expected 8{H*) to sustain 
its low temperature value up to the temperature where 
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1 Tesla 

■ fit with T S =37.3K±1.7K 

2 Tesla 

• fit with T S = 21.9K±2.5K 

3 Tesla 

■ fit with T c - 14.4K±1.0K 



FIG. 5: The decay of 5 as a function of the annealing tem- 
perature T a . The values of (S up (T a ) — Sf.c.) are normalized 
by the value of (8 U p(lAK) — Sf.c.)- Solid Lines are fit to 
exponential decay of first order. With T c = 89K The fitted 
order of decay fs. jg 0.16 ± 0.01, 0.25 ± 0.03 and 0.42 ± 0.02 
for 1, 2and3 Tesla respectively. 





Bean-Livingston 




Nodal surface 


* 


T s 


o 






T/T 



FIG. 6: Field of first vortex entry as a function of the reduced 
temperature ^- . The lines correspond to the Bean-Livingston 
field of first vortex entry, using H c = 1 Tesla [12], in the case 
of conventional superconductor (dashed line) and in the case 
of nodal surface calculated by Iniotakis et. al. @] (full line). 
The circles are represent the field where the hysteresis value of 
the splitting, (S up (T a , H*) — Sf.c.) is at maximum that is Hm 
taken from Fig. [4] The lowest temperature (noted also with a 
question mark) is only a lower boundary since no maximum is 
observed. The stars represent the decay order T S {H*) taken 
from Fig. [5] Both methods of analysis show the enhancement 
of H s for low temperatures. 



H* is equal to the vortex entry field. 

There are several possible reasons for the observed 
broad maximum (Fig. [4]) and for the continuous decay of 
S up (Fig. [5|). A first one is the surface roughness of the 
sample that can have weak spots for vortex entry where 
the surface is not (110) oriented. A second one is the 
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height of the barrier which may become of order KbT 
at high fields. One must also take into account the finite 
thickness of the sample which is only a few penetration 
depths. 

The formation of an id xy subdominant order parame- 



ter 



ToL 13l | , can be expected to reduce the paramagnetic 
currents and hence H s . However, at fields of less than 
3 Tesla the id xy component disappears above about 7K 
[l3j ]. Therefore, it can not effect the results presented in 
Fig. [4] and Fig. [5] at temperatures above 14K. As for the 
1.4K data (Fig. QJ, it only allows to determine a mini- 
mum value for H s which was not calculated by Iniotakis 
et.al. 0. 

To summarize the results, we show Fig [6] values of 
— - and of HM{qT-)- Within error bars, both sets 
of values fall on the same line, showing that both cap- 
ture the same physics. The dashed line in Fig. [6] is 
the temperature dependence of the Bean-Livingston field 
of first entry in a conventional superconductor, taking 
H c = 1 Tesla pjj]. The solid line is the temperature de- 
pendence of the Bean-Livingston field of first entry in the 
case of a nodal surface, taken from Iniotakis et. al. 
again with H c = 1 Tesla. The predicted large enhance- 
ment is clearly seen experimentally at low temperatures. 

The reason for the low temperature enhancement of 
the field of first entry is easily understood. As shown by 
Fogelstrom et.al. Q the Andreev - Saint-James states 
carry paramagnetic currents at nodal surfaces. Iniotakis 
et.al. 3] have shown that these currents strongly increase 
at low temperatures, thus reducing the usual repulsive 
Lorenz force due to Meissner currents. This reduction 
of the repulsive force allows the attractive force of the 
image vortex to remain dominant up to higher fields at 
low temperatures. 

In order to better understand the barrier further the- 
oretical and experimental work is needed in lower tem- 
peratures regime where the barrier is strongly enhanced. 
The doping dependence of the enhancement and espe- 
cially the appearance of a subdominant nodal order pa- 



rameter can strongly effect the Bean-Livingston barrier 
enhancement due to the node removal effect. In such a 
case, strong over-doped samples and/or high magnetic 
fields are necessary in order to allow the existence of the 
subdominant order paramter at temperatures as high as 

0.3T C [m m. 
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